****** Table_11_1: Specification tests
use data/Demo_book_country_3.dta, clear
************ democracy variables

* polyarchy, 100-point scale
gen v2x_polyarchy_100 = (v2x_polyarchy / .924) * 100

* polyarchy, missing obs imputed
impute v2x_polyarchy lexical_index inc_chall_combined , gen(v2x_polyarchy_imp)
* 100-point scale
replace v2x_polyarchy_imp = v2x_polyarchy_imp * 100
gen v2x_polyarchy_imp_100 = (v2x_polyarchy_imp / 92.4) * 100
replace v2x_polyarchy_imp_100 = . if sovereign_Cojocaru ==.


*1) Bivariate
qui reg v2x_polyarchy_imp_100 eur_pct_est_smooth i.year, vce(cluster country_id)
qui gen sample = e(sample)
qui tab country_id if sample == 1, nofreq
qui local countries=r(r)
qui tab year if sample == 1, nofreq
qui local years=r(r)
qui outreg2 using output/Table_11_1.doc, nor noni nonot nocon stats(coef tstat) addstat("Countries", `countries', "Years", `years', "R-squared", e(r2)) addn("*** p<0.01, ** p<0.05, * p<0.1") se bdec(3) keep(eur_pct_est_smooth) replace
qui drop sample

*2) Geography
qui reg v2x_polyarchy_imp_100 eur_pct_est_smooth predallports53_100km lp_lat_abst i.year, vce(cluster country_id)
qui gen sample = e(sample)
qui tab country_id if sample == 1, nofreq
qui local countries=r(r)
qui tab year if sample == 1, nofreq
qui local years=r(r)
qui outreg2 using output/Table_11_1.doc, nor noni nonot nocon stats(coef tstat) addstat("Countries", `countries', "Years", `years', "R-squared", e(r2)) addn("*** p<0.01, ** p<0.05, * p<0.1") se bdec(3) keep(eur_pct_est_smooth) append
qui drop sample

*3) Additional geographic covariates
qui reg v2x_polyarchy_imp_100 eur_pct_est_smooth predallports53_100km lp_lat_abst tropical avg_temp frstdays Island landlock_new rivdist_abs_new rugged mountains c.prec_gpcc##c.prec_gpcc irri_impact5 soil suita desert i.year, vce(cluster country_id)
qui gen sample = e(sample)
qui tab country_id if sample == 1, nofreq
qui local countries=r(r)
qui tab year if sample == 1, nofreq
qui local years=r(r)
qui outreg2 using output/Table_11_1.doc, nor noni nonot nocon stats(coef tstat) addstat("Countries", `countries', "Years", `years', "R-squared", e(r2)) addn("*** p<0.01, ** p<0.05, * p<0.1") se bdec(3) keep(eur_pct_est_smooth) append
qui drop sample

*4) Kitchen sink specification
qui reg v2x_polyarchy_imp_100 eur_pct_est_smooth predallports53_100km lp_lat_abst murdock_dem chrstprotpct Muslim durbritain_new indep_yrs_since_ln statehiste1500_01n al_ethnic al_language al_religion e_Total_Oil_Income_PC_ext lpd1500 el_settlermort el_ME agritran i.e_regionpol4 i.year, vce(cluster country_id)
qui gen sample = e(sample)
qui tab country_id if sample == 1, nofreq
qui local countries=r(r)
qui tab year if sample == 1, nofreq
qui local years=r(r)
qui outreg2 using output/Table_11_1.doc, nor noni nonot nocon stats(coef tstat) addstat("Countries", `countries', "Years", `years', "R-squared", e(r2)) addn("*** p<0.01, ** p<0.05, * p<0.1") se bdec(3) keep(eur_pct_est_smooth) append
qui drop sample

*5) Full specification: including all factors that seem - on theoretical and empirical grounds - to below
qui reg v2x_polyarchy_imp_100 eur_pct_est_smooth predallports53_100km lp_lat_abst murdock_dem chrstprotpct Muslim durbritain_new i.year, vce(cluster country_id)
qui gen sample = e(sample)
qui tab country_id if sample == 1, nofreq
qui local countries=r(r)
qui tab year if sample == 1, nofreq
qui local years=r(r)
qui outreg2 using output/Table_11_1.doc, nor noni nonot nocon stats(coef tstat) addstat("Countries", `countries', "Years", `years', "R-squared", e(r2)) addn("*** p<0.01, ** p<0.05, * p<0.1") se bdec(3) keep(eur_pct_est_smooth) append
qui drop sample
